Wigner crystallization of photons in cold Rydberg ensembles 
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The coupling of weak light fields to Rydberg states of atoms under conditions of electromagnetically induced 
transparency (EIT) leads to the formation of Rydberg polaritons which are quasi-particles with tunable effective 
mass and long-range interactions. Confined to one spatial dimension their low energy physics is that of a 
moving-frame Luttinger liquid which due to the long-range character of the repulsive interaction can form a 
Wigner crystal. We calculate the Luttinger K parameter using density-matrix renormalization group (DMRG) 
simulations and find that under typical slow-light conditions kinetic energy contributions are too strong for 
crystal formation. However, adiabatically increasing the polariton mass by turning a light pulse into stationary 
spin excitations allows to generate true crystalline order over a finite length. The dynamics of this process and 
asymptotic correlations are analyzed in terms of a time-dependent Luttinger theory. 

PACS numbers: 32.80.Rm,42.50.Gy,32,80.Qk 



The extraordinary properties of Rydberg atoms Q), such as 
large dipole-dipole interactions and long life-times, are cur- 
rently attracting much attention. The interest ranges from 
quantum information ||2}5) to many-body phenomena 15 HT31 . 
So far only few works considered the effect of interactions 
onto the light fields ifTBT - EOl . In recent experiments [16] it 
was shown that under EIT conditions the Rydberg interaction 
leads to a non-local, and strongly non-linear behaviour of the 
probe field [22 23|. This gives rise to, e.g., the formation of 
a small avoided volume which contains at most one excitation 
|20| . which can be described by considering only two excita- 
tions [3|. In the present paper we want to explore the many- 
body properties on larger length scales. One of the simplest 
but most dramatic effects resulting from a long-range repul- 
sive interaction is the formation of a Wigner crystal, predicted 
for electrons in the early days of quantum mechanics [25 1. We 
will show that a similar phenomenon can be observed in a di- 
lute ID gas of photons coupled to Rydberg atoms, which is not 
possible in e.g. Kerr-type point-interacting systems I26l l27l 



Under conditions of EIT and small excitation densities, the 
coupling between photons and Rydberg atoms leads to the for- 
mation of dark-state polaritons (DSP) [28, 29], which behave 
like massive particles with strong repulsive and long-range 
interaction. We analyze the formation of a quasi-crystalline 
state of polaritons in ID using DMRG simulations and time- 
dependent Luttinger-liquid (LL) theory. We show that under 
slow-light conditions the moving-frame ground-state displays 
density-wave correlations which decay very fast in propaga- 
tion direction due to the small polariton mass. However, con- 
verting them into stationary spin excitations by decelerating a 
light pulse to a full stop inside a gas of Rydberg atoms II301I3T1 
can lead to a state with perfect crystalline order over the length 
of the medium. 

To be specific, we consider an ensemble of N atoms with 
a three-level linkage-pattern [cf. Fig. [IJb)], composed of 



a ground-state |g), intermediate state |e) and metastable 
Rydberg-state |r). The transition |g)— |e) is driven by a quan- 
tized probe field E = J^£(r, t) e - i ( w p*- c 'p r ) + H.a., with 

carrier frequency lj p and wave-vector q p . £(£') are normal- 
ized field amplitudes corresponding to annihilation (creation) 
of a photon and are slowly varying in space and time. The 
transition |e) — |r) is coupled via a strong external control field 
with Rabi frequency il, carrier frequency lo c and wave-vector 
q c . We chose the z-axis as the common propagation direction 
of the probe and control field and define the one- and two- 
photon detunings as A = u c — cj g — uip, 5 = u> r —u g —oj c —ujp, 
respectively, where w gcr are the energies of the correspond- 
ing atomic states (h = 1). 

In the absence of Rydberg interactions the system Hamilto- 
nian can be diagonalized using adiabatic eigen-solutions, the 
dark- and bright-state polaritons (BSP), which fulfill approxi- 
mate bosonic commutation relations ll28l|29l . Following ll32ll 
we define the DSPs as ^ = cos 0£ 
BSPs as <l = smO£ 



sm9E gr e l( -^+^ z ,and 



cos6»E gr e^+«p) 2 where tan 2 9 = 
g 2 n/n 2 . Here E M „ are continuous atomic spin flip opera- 
tors, and n is the atom density. The DSP propagates loss- 
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Figure 1 . (a) Schematic setup for the creation of dark-state polaritons 
in a medium on length L. (b) Effective atomic linkage pattern for EIT 
in Rydberg gases. The weak quantized field £ is off-resonantly driv- 
ing the |g) — |e) transition with a one-photon detuning A, whereas 
the strong control field fi is driving the |c) — |r) transition with a 
final two-photon detuning S. 



less with group velocity v g — c cos 2 9, while the BSP has 
a velocity c sin 2 9 and is subject to losses with rate Sl 2 /r, 
where T = 7 + i A, with 2j being the spontaneous decay rate 
of state |c), and il 2 e = g 2 n + il 2 . Near single-photon reso- 
nance, i.e., |A| < 7, and for an optically thick medium, i.e., 
L 3> £ a bs, where i a b s = g 2 n/c-f is the resonant absorp- 
tion length in absence of EIT, an input bright-polariton will 
quickly be damped out. In the following we will consider the 
opposite case |A| >• 7, where absorption is irrelevant. How- 
ever if cos 2 9 <C sin 2 9 w 1 and if light pulses of finite length 
are considered an input bright-polariton can be disregarded as 
it will quickly escape the medium. Thus the BSP can be elim- 
inated and the free dynamics of the DSPs is governed after a 
short transient by [ 33 J 
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where p z = — id z , p^ = — iVj_, and <5(r) = 8(r) sin 2 9. This 
corresponds to an effective Schrodinger equation for particles 
with tensorial mass and additional drift term, moving in an 
external potential <5(r) sin 2 9. The drift is determined by the 
EIT group-velocity v g , and the masses are mj 1 = v g -L a bs 
and mj 1 = ^- fi"l l32l . The above model is valid as long as 
the BSP amplitude is negligible and [33] 
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The first condition describes the regime of perturbative cou- 
pling between DSP and BSP, the second denotes the region of 
slow-light dispersion lfTl l29ll . L c har is a characteristic length 
scale of the DSP. 

Let us now take into account interactions between the atoms 
in their Rydberg-state |r), H Ryd = |£\ .._., ai- 3 / V dd (r 3 - 

T , j )a^' ) , where V dd (r) = C a /\r\ a with a = 6 (a = 3) being 
the van-der-Waals (resonant dipole) interaction. In the con- 
tinuum limit and transforming to polaritons, we find to lowest 
order in cos 9 



The effect of the interaction is equivalent to a two-photon de- 
tuning. Consequently, the interaction shift has to be smaller 
than tt 2 /\ A I which can be translated into a minimal distance 
a c = (C Q 7/ri 2 ) 1 / Q the DSPs have to keep to ensure the va- 
lidity of the model. As shown in [3| for the case of a res- 
onant interaction (i.e., A = 0) and large optical depth, an 
incoming coherent light pulse will quickly develop strong 
anti-bunching with a minimum separation length in propa- 
gation direction corresponding to the EIT blockade radius 
0>h = (Caj/Cl 2 ) 1 '". A similar effect happens for A ^ 
due to the fast escape of the BSP. Since under slow-light con- 
ditions cos 9 = £l/n e <C 1 the initial preparation produces 
DSPs with a mutual distance larger than the critical value a c 
and a vacuum of BSPs. Consequently the system can be well 
described by H = Hq + H mt . 



To address the question whether the interaction leads to 
Wigner-crystallization of polaritons we restrict ourselves to 
one dimension (ID). This can be achieved e.g., by using elon- 
gated cigar-shaped atomic ensemble with transverse extent 
smaller than the blockade radius [20 1, or atoms in hollow-core 
fibers [35 36 1 or trapped in the evanescent field of ultra-thin 
optical fibers ll37l [38l . The low-energy physics can be de- 
scribed in terms of a Luttinger liquid (LL) [ 39 1 . The LL model 
allows for an exact treatment also in the case of bosons Il40l 
with l/\x\ "-interactions, as long as a > 1, Transforming to 
a frame co-moving with the EIT group-velocity removes the 
drift term, oc v g d z , in eq. ([TJ. Assuming a fixed number of ex- 
citations and working off single-photon resonance | A | 3> 7, 
we follow the standard LL-approach [39| to construct an ef- 
fective low-energy Hamiltonian 
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IT and are conjugate fields with <f>(x),il(y) = iS(x — y). 
u and K are the sound velocity and the Luttinger parame- 
ter, respectively. The X-parameter governs the long-range 
behavior of the charge-density-wave correlations (CDW) in 
the ground-state. E.g., the oscillatory part of the density cor- 
relations is given by (p(z)p(0)) OBC ~ p 2 , cos(2n p z) z~ 2K . 
Here p(z) = W{z)^{z) and p is the average density. As 
first-order correlations decay as (&(z)ty(0)) ~ Z ~ X ^ 2K ^ the 
point K = 1/2 marks the crossover from a regime where 
superfluid order dominates (K > 1/2) to a regime with pre- 
dominant CDW correlations of period 1/po (K < 1/2). 

At this point it should be noted that one can create true 
crystalline order by adding a weak periodic lattice poten- 
tial S(x) = Sq sm(2irx/d), which leads to a sine-Gordon 
Hamiltonian [39| for commensurate fillings p = l/(sd), 
s € N. This model exhibits a quantum phase transition to 
a gapped ordered phase for arbitrarily small but finite 6q, if 
K < K s = 2/s 2 J39j ED. To avoid the necessity of a co- 
moving lattice potential one then should consider stationary- 
light polaritons lHI l42l l43l . 

Although no exact expression for K exists, an approximate 
closed formula was given in [44|: 
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Here C,(x) is the Riemann Zeta-function and we introduced 
the dimensionless interaction strength 9 = Po ~ 2 ™^° ■ 

To check this expression we determined K numerically us- 
ing DMRG [45 1 and made use of the fact that K/u = ^p\x is 
determined by the compressibility x^ 1 — Po^~ — PqL^t 
ll39l . Furthermore uK — Trpo/m, which is true for any 
Galilean-invariant model [46] . We have validated the nu- 
merical procedure for the case of the integrable Lieb-Liniger 
model, where \ can t> e calculated exactly as a function of in- 
teraction strength. Using a proper discretization of the model 
B71 leads to the results for K shown in Fig. [2] One recognizes 
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Figure 2. Luttinger parameter K as a function of interaction strength. 
The continuous (dashed) line shows the analytical approximation 
1441 for a — 6 (a = 3). Diamonds (squares) show results from 
DMRG calculations with open {red} and periodic {blue} bound- 
ary conditions (BC). The numerical parameters were Ax — l/10po 
with (i) L — 20/po and bond dimension d=16 for open BC or (ii) 
L — 10/ po and d=32 for periodic BC. Quantum Monte Carlo results 
for a = 3 can be found in 1401 . 



good agreement for dipolar interactions {a = 3) but a shift 
towards smaller values of K for vdW interactions (a = 6). 
Nevertheless we will use the analytic expression in the fol- 
lowing. 

Let us first discuss the time -independent case D,(t) = O = 
const. We concentrate on the ground-state where CDW corre- 
lations should be most pronounced. In Fig. [3]we have plotted 
the normalized two-particle correlation g^Jz) in the ground 
state of H = H + H; nt obtained by DMRG corresponding to 
different values of K. The large-z behavior is well described 
by LL correlations and one recognizes well pronounced os- 
cillations for K < 1/2. For small distances the plots show 
an extended spatial region over which g( 2 )(z) vanishes. The 
correlations around z = become more suppressed than in 
the case of free fermions, which is the strongest possible for 
point interactions [26 48 1. Thus the long-range interaction, 
eq. (J3), prevents that two DSPs approach each other which 
also proves the consistency of the effective model. 

We can use eq. ([5} to estimate the critical interaction 
strength required to enter the CDW dominated regime, i.e., 
K < 1/2, giving r/ a ® a -i t = 3. 9 is proportional to the effec- 

2, 

tive mass of the polaritons m ~ v~ ~ ^r which is different 
along longitudinal (mil) and transverse directions (mi), and 
can be tuned via the control field Vt, Considering vdW inter- 
actions (and setting rj e = 1) we find for m = m\\ . 
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where we have introduced the optical depth per (initial) block- 
ade radius OD^ = ab/L a t>s- Note that condition (|2| translates 
into p L ahs < 7/|A|. Using e.g. p L ahs = j/\A\ = 1/100, 
we find that the optical depth per blockade at cr i t has to be 
OD^ > 30. This shows that a crystalline structure will be 
challenging to prepare along the propagation direction. Along 
the transvere direction the conditions are more relaxed. To 



overcome this difficulty we will now discuss an alternative 
based on light storage. 

Let us consider an initial polariton pulse close to the ground 
state in the co-moving frame and time-dependent control 
fields. Decelerating the DSPs by reducing v g in time pre- 
serves their spatial structure (and density p ), in the absence 
of interactions ||29l . At the same time their effective mass is 
increased, which suggests a decreasing K according to eq.|5]) 
for interacting DSPs. When the pulse is brought to a complete 
stop, K(t) approaches zero potentially leading to true long 
range order. If v g is switched to zero instantaneously, the ini- 
tial spatial correlations will be frozen. Thus the switching has 
to be done smoothly on a time scale r long enough for corre- 
lations to propagate through the system. The latter process is 
determined by the speed of sound u(t) — irpo/(m(t)K(t)). 
For small K we find the scaling K ~ 1/ \fm ~ ^/Vg, i-6., 
the speed of sound decreases only with the square root of the 
group velocity, u(i) 



In order to describe the adiabatic switch-off we consider 
the LL Hamiltonian (p| with time-dependent parameters K (t) 
and u(t) [49|. Using the standard decomposition of the fields 
into bosonic momentum modes b p , M, ||39l yields 



1 TTPn 
H = - 
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\p\{w(t)blb p -?f(blbl p + b- p b p )}, 



2 m(t) -V 

(7) 
= 1 + 1/K(t) 2 , g(t) = 1 - l/K{tf and 
V6C'6Po rn (t)/'^' K ) 1 ^ 2 ^ according to eq. |5|. 



where w(t) 
K(t) = (1 

The corresponding Heisenberg equations of motion can be di- 
agonalized by a time-dependent Bogoliubov transformation. 
Since the transformation matrix is itself time-dependent this 
leads in general to an infinite hierarchy. However, if K(t) 
and u(t) are chosen such that K(t) ex K(t) u(t) the hierar- 
chy truncates after the second iteration. This requires m(t) = 
mof(t) = m e x W sinh [x(t)] e - arcosh ( c ) / \/ C 2 - 1, with 
x(t) = arcosh(i/r + C*), where C = (K$ + 1)/(2K ), 
K Q = K(t = 0) and r being the characteristic deceleration 
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Figure 3. Normalized two-particle correlation, main panel: 1 — 
g( 2 '(z), in double logarithmic scale, inset: g^ 2 '(z) in linear scale. 
Full lines show numerical results for interaction strength increasing 
from orange over black to blue. The particles are subject to periodic 
boundary conditions. Note that at poz > 5 finite size effects become 
noticable. 



time. The time-dependence of the control-field can then be 

computed as O(i) = gy/n[f(t)c/v g (Q) — l) , where 
too = m(t = 0) and m(t)~ l = v ff (i)L a b s A/7 was used. Then 
K(t) — e~ x ( t ) and we find the oscillatory part of the equal 
time density-density correlations given by 
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where £(t) = y/lftp 2 - Hn(K(i)/K ). The first term ~ 
z -iK(t) (j escr ibes the adiabatic behavior, where the spatial 
dependence of correlations is determined by the instantaneous 
value of K(t). The exponential describes non-adiabatic cor- 
rections. As a consequence of the finite speed of sound, these 
corrections lead to a "crossover" from the power law with adi- 
abatic exponent K(t) to one with the initial exponent Kq at 
the characteristic distance Iq = r/^CepQ r, as can be seen from 
Fig. HI Maximizing the switch-off time r results in large spa- 
tial regions lo with quasi-crystalline order. On the other hand, 
r has to be sufficiently small to bring the pulse to a complete 
stop after a length L^ = L dtv g (t) less than the medium 
length L. Setting Ld = L and using the above protocol, we 
find 
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Noting that Kq is close to unity and that poLahs < 7/| A|, this 
fraction can approach unity showing that a crystalline order 
over the whole medium is possible. 

It should be noted that changing the control-field in time 
leads to an additional coupling between the DSP and BSP 
(29). The decay rate due to this coupling is given by jg = 
r y8 2 /g 2 n and we need J dtjg(t) <C 1. Using the above pro- 
tocol we find cr/L. dbs > 4K$/(K$-1) 2 . For K « 0.99 and 
iabs ~ 5 pm we have r 3> 40 ps, which is certainly feasible. 

So far we have assumed that the initial state for the light 
storage is the moving-frame groundstate of the Hamiltonian 
Q. Let us now discuss the effects of initial excitations. As 
the system is non-integrable it is reasonable to assume that the 
state of the DSPs after the initial preparation is thermal (we set 
kg = 1). In a thermal state all correlations decay exponen- 
tially for distances larger than Lt ~ u/T, with correlation 
length £ T - L T jK = irp /{mTK 2 ) ED. Since for a non- 
interacting gas and adiabatic mass changes T(t) ~ l/m(t) 
holds, one naively expects that the correlation length increases 
£ T rv 1/K 2 (t) as long as one stays adiabatic. Evaluat- 
ing the bosonic correlation functions in eq. dHl with a ther- 
mal distribution we find a slightly different result (cf. Fig. 
p). Correlations decay as exp(— \z\ 2 / L 2 orT ) for length scales 
\Z\ < icorr = 2^/lol^/7rKo x {\u{K / 'K{t))) x l A which 
crosses over to an exponential decay for larger length scales. 
Here L^ is the initial thermal correlation length. 

To estimate the maximal initial temperature of the polariton 
system we observe that any polariton components with fre- 
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Figure 4. main panel: space dependent amplitude of the oscilla- 
tory part of the correlation function in the long time limit, Kg = 
0.8, K (t) — 5 x 10 , The dashed blue line shows the spatial decay 
of density-density correlations for zero Temperature, which shows a 
crossover from adiabatic to diabatic algebraic decay at lopo = 100 
(indicated by rightmost vertical line). The solid red line shows the 
modified decay for initial Temperature corresponding to a thermal 
length Lrpo = 10 (leftmost vertical line) which shows a crossover 
to exponential decay at length scale L COIX po ~ 40. 
inset: Q(t)/Q(0) for K = 0.8 and v gT (0)/c = 10" 5 . 



quencies larger than the off -resonant EIT line-width f2 2 /|A| 
will escape [33 1. Thus we can give as a reasonable estimate 
for an upper limit T 



hjM- With this we find 






PoL- d 




(10) 



which gives a measure for the final correlation length in units 
of the total length. Although the first term on the right side is 
less than unity, the whole expression can approach unity. 

In summary we showed that the combination of EIT with 
interacting Rydberg gases leads to strongly interacting light- 
matter particles, termed Rydberg polaritons. We discussed 
the experimental requirements needed to obtain a quasi-long- 
range-ordered ground-state corresponding to a moving-frame 
Wigner crystal of Rydberg excitations in ID by mapping the 
problem to a Luttinger liquid. Numerical and analytic re- 
sults showed that under slow-light conditions the kinetic en- 
ergy contributions in the longitudinal direction are too large 
to enter the density-wave dominated regime. Using a time- 
dependent Luttinger liquid approach we showed, however, 
that decelerating a light pulse in a gas of Rydberg atoms to 
a full stop over a sufficiently long deceleration time can cre- 
ate true crystalline order over a substantial fraction of the 
medium. Turning the Wigner crystal of spin excitations back 
into electromagnetic fields by a sudden switch-on of the drive 
field l28l produces a train of photons with long range crys- 
talline order. 
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WIGNER CRYSTALLIZATION OF PHOTONS IN COLD 
RYDBERG ENSEMBLE - SUPPLEMENTAL MATERIAL 

In this supplementary we derive the effective Hamiltonian 
(1) of the main text in one spatial dimension and explicitly 
establish its conditions of validity. The treatment is based on 
the formalism developed in [ 1 1. The interaction of probe and 
control fields with the three-level atoms shown in Fig. 1 of 
the main text in the absence of Rydberg interactions can be 
described by the following atom-light coupling Hamiltonian 
in a rotating frame 

// = /d 3 r{ A± ee (r) + 6± rr (r) (11) 

- £ffi re e lq <= r + 5 V^(r)S eg (r)e iq " r + H.a.}, 



where all quantities are defined as in the main text and 
(v\/\fA~N are continuous atomic flip 



E M „(r 



Z^jeAV \f i )jj\ 
operators defined on a small volume AV(r) centered around 

position r containing AiV 3> 1 atoms. Assuming that all 
atoms are initially prepared in the ground state \g) and con- 
sidering weak probe fields, i.e. a photon density much less 
than the atom density, we can treat the light-atom coupling 
perturbatively. Consequently, in lowest order of the atom-field 
coupling g we find the Heisenberg-Langevin equations for the 
atomic operators 



at ge ~ 



at 9r 



- (*A + i)t ge + F ge 

+ igy^£e l ^ r + irr£ ffr e-' iq <= r , (12) 

- i5tg r + iVtt ge e^ T . (13) 



Here the F ge is a delta-correlated Langevin noise operator 
associated with the decay from the intermediate (excited) 
state |e) which is necessary to preserve commutation rela- 
tions |2|. One easily verifies that the correlation functions 
of the Langevin operators are proportional to the population 
in the excited state |e). In the linear response and for suffi- 
ciently small two-photon detuning 8 this population is small 
and we can safely ignore the noise operators in the following. 
If need be these operators can be re-introduced by hand using 
the fluctuation-dissipation theorem. To arrive at a closed de- 
scription of the atom-field system we also need the equation 
of motion for the slowly varying probe-field envelope £ (r, t). 
Restricting ourselves to a one-dimensional problem the dy- 
namics of the probe field is described by a truncated wave- 
equation in paraxial approximation 



d_ 

at 



d_ 

dz 



£(z,t) = ig^t ge {z,t)e^ T 



(14) 



Transforming Eqs. ( fT2"p4} into Fourier space according to 

f(z, t) — j dk e~ tkz f(k, t) yields the following matrix equa- 
tions 



dt 



X 



-MX 



(15) 



where X T = {£, t gr e i ^+^ T , S 3e e' q p r } and the Hamilto- 
nian matrix reads 



-fee 





~9\/n 





5 


-n 


g\pn 


-Q* 


A — 17 



Changing the basis to a description in terms of dark- and 
bright-polaritons Y T = {*, $, Egge" 1 " 1 "} via * = cos8£ - 



sm9T, gr e l( -^+^)r and $ 



sin9£ + cos6>£ 9r e i ( c 'p+ c ^ r 



yields the equation of motion d t Y 
H' = 



5 sin 2 6- kc cos 2 9 -sin9 cos 9(6 + kc) 
— sin 9 cos 9(5 + kc) 5 cos 2 9 — kc sin 2 9 —fl e 
-O e A - 17 



Assuming that the time evolution is slow compared to | A — ij\ 
we can adiabatically eliminate the optical polarization E ffe 
which yields the coupled equations for bright and dark po- 
laritons 



d 
dt 



= -m' 



(16) 



with 



H" = 



S sin 2 9 — kc cos 2 9 
— sin 9 cos 9 (6 + kc) S cos 2 9 — kc sin 



— sin 9 cos 9(8 + kc) 

A — 17 



For A > and under slow-light conditions, i.e. sin 9 3> 
cos 2 9, one recognizes that the off-diagonal coupling terms 
are small compared to the difference of the diagonal elements. 
Under these conditions the dynamics of dark and bright po- 
laritons approximately separates and one can treat their cross 
coupling perturbatively. Within this perturbative treatment the 
effective equation of motion of the dark polariton *£> up to sec- 
ond order of the off-diagonal coupling is given by 



dt 



* = -i (S sin 2 9 - kc 



ccos 



* 



sin 2 6>cos 2 6»((5 + fcc) 



(5 + kc) (sin 9 - cos 2 ( 



Q2 

A — ?7 



-*. 



To arrive at an even more simplified but more transparent 
equation we assume 6 > and require 5 + kc -c fi 2 /|A| 
for all relevant values of fe, which implies in particular 



0<<5< 



IAI 



(17) 



where we used | A| ^> 7, In this limit we find that the dynam- 
ics of the DSPs is described by 



dt 



-* = -id 1 



8 A cos 2 9 



02 



* + ikv„ 



-'£'* 



v„cA 



fc 2 *, 



(18) 



where we approximated sin 9 « 1, The first term on the 
right hand side describes an energy offset due to a finite two- 
photon detuning. The second term accounts for the propaga- 
tion with group velocity v g = c cos 2 9. The third term de- 
scribes the quadratic dispersion with effective mass mj 1 = 
c 2 cos 2 #|A|/f> 2 «t> g L abs |A|/ 7 . 

Condition ([TTJi also determines the validity of the interac- 
tion Hamiltonian (3) of the main text. As already pointed 
out in the main text, the Rydberg interactions effectively in- 
duce a space-dependent two-photon detuning [ 3 1 which com- 
bined with eq. ( fTTj ) leads to the critical minimal distance 

a c = {c al /ni)y^ 

We can also interpret the validity of the perturbation theory 
as a condition of a maximal fc-value until which a separation 
into DSPs with slow-light dispersion and fast moving BSPs is 
valid. This condition reads 



\kc\ < 



ii. 



(19) 



leading to |/e max | = £7 2 /|A|c. Plugging this into the zeroth 



order dispersion relation of the DSP leads to the maximal en- 
ergy Wmax = w g fc max = £7 2 /| A|, which in the end determines 
the maximal temperature of the DSP gas. 

Finally estimating the typical A; -value of the system via the 
inverse characteristic length scale, i.e. k ~ 1/L c h ar we can 
rewrite condition (fT9ll and obtain 



L 



ab.s 



< 



J char 



7 
IAI' 



(20) 



which is just condition (2) of the main text with the appproxi- 
mation sin (9 ss 1. 
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